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1. INTRODUCTION 

In the Standard Model (SM), the flavor-changing neutral processes are very rare and are of 
obvious interest in the search for new physics. In the Flavor-Changing Neutral Current (FCNC) 
decays D° — > £ + £~, where £ is an either electron or muon, are strongly suppressed by the Glashow- 
Iliopoulos-Maiani (GIM) mechanism |JI]]. These decays cannot occur at tree level in the SM. The 
branching fraction of the decays D° — > £ + £~ are predicted to be i^(10~ 13 ) [gj. 

Most of the attention on FCNC decays has been focused in the K and B meson sectors and 
less in the charm meson sector. It is due to the fact that the SM expectations for D° — D° mixing 
are very small compared to the K° — K° and B° — B° mixing. However, the FCNC decay in the 
charm sector is unique due to its decays involve an up-type quark which implies into an effective 
GIM cancellations and new physics. The decay modes of Lepton-Flavor Violating (LFV) which 
corresponding to two leptons with two oppositely charged of different flavor and Lepton-Number 
Violating (LNV) decays where two leptons have the same charge are forbidden in the SM. 

Figure [j] shows the Standard Model short-distance contributions to the c — > u£ + £~ transition. 
The branching fraction for the decay of D — > X u i + t~ is predicted to be ^(10~ 8 ) [g, [3J]. The decay 
of c — > u£ + £~ is screened by the long distance contributions. It is also expected to dominate over 
the short distance contributions in D° — D° mixing. The long distance contributions were shown 
to be largely dominant in c — > u£ + £~ . The experimental upper bounds on the branching fraction 
of c — > u£ + £~ is presently in the range of <??(10~ 5 ) [Q]. It is an order of magnitude larger than 
the Standard Model prediction for specific channels [Bh. The highest rate of D — > V£ + £~ channel 
with V = p, (0, (j),K* is the decay of D+ — > p + £ + £~. It is predicted at the highest rate m 3 x 10~ 5 , 
however the there are unfortunately no experimental data on this channel. 




Figure 1: Standard Model short-distance contributions to the c — > u£ + £ transition. 



2. The BABAR DETECTOR AND DATA SET 

The BABAR detector was operated at the PEP-II asymmetric-energy storage rings at the SLAC 
National Accelerator Laboratory. The data used in this analysis were collected with the BABAR 
detector. The data sample consist of an integrated luminosity of 384 fb~ x for X c + — > h ± £ T £^ + 
and 468 fb^ 1 for D° £+£~ accumulated at the T(4S) resonance and 40 MeV below the T(4S) 
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resonance. The asymmetric energy of the PEP-II e + and e beams result in a Lorentz boost /3y « 
0.55 of the BB pairs. 

A detail description of the BABAR detector is presented elsewhere [§]. The momenta of the 
charged particles are measured in a tracking system consisting of a 5-layer double sided silicon 
vertex tracker (SVT) and a 40-layer drift chamber (DCH). The SVT and DCH operate within a 
1.5 T solenoid field and have a combined solid angle coverage in the center of mass frame of 
90.5%. A detector of internally reflected Cerenkov radiation (DIRC) is used for charged particle 
identifications of pions, kaons, and protons with likelihood ratios calculated from dE/dx measure- 
ments in the SVT and DCH. Photons and long-lived neutral hadrons are detected and their energies 
are measured in a CsI(Tl) electromagnetic calorimeter (EMC). For electrons, energy lost due to 
bremsstrahlung is recovered from deposits in the EMC. 

3. ANALYSIS 

We select charm hadron candidates X c with center of mass frame momentum greater than 2.5 
GeV/c to suppress combinatoric background. The e + e~ invariant mass is required to be greater 
than 200 Mev/c 2 in order to reject photon conversion and 71° decays to e + e~y. For the Dt% — > n(j>, 
<j> — > £ + £~ decay mode, we excluded events with 0.95 < ni(e + e~) < 1.05 GeV/c 2 and 0.99 < 
m(n + iJ.~) < 1.05 GeV/c 2 to reject the decays through the resonance. The QED backgrounds 
was suppressed by requiring at least five tracks in the event and that hadron candidate be consistent 
with the electron hypothesis. After the initial event selection, significant combinatorial background 
contribution, we use three discriminating variables in likelihood ratio: charm hadron candidate, 
total reconstructed energy in the event, and flight length significance. 

To measure the signal events we use extended, unbinned, maximum-likelihood. These signals 
are converted to the known charm branching fractions by normalization. To reduce the systematic 
effects we choose normalization modes with kinematics similar to the kinematic of the signal de- 
cays. For decays of D + and D+ mesons, the normalization mode is n + <j) where — > K + K . For the 
decays of A+, we choose the decays of A+ — > pK~n + . Figures ||-|| show the fitting results of the 
invariant mass of X c + — > h ± £ T £^ + decays. The dashed curves show the background components for 
the dimuon modes in which muon candidates arise from misidentified hadrons. Detail information 
on the likelihood selection, fitting procedure, systematic uncertainties, and fit results are available 
here [0]. We calculate the upper limits on the ratio of the branching factions at 90% confidence 
level (CL): Sipf. -> n ± £^£^ + )/^(D+ ) -> K + (j>), 38(0+. -> K ± e^ + )/^{D+ -> 71+0), and 
J?(A+ — > p(~)l^ffl )/38(A+ — > pK~n + ). The most significant signal is seen in the decay of 
A+ — > PH + IJ.~ with yield of 11.1 ±5.0(stat) ±2.5(syst). It has a statistical significant of 2.6(7. It 
is corresponding to 90% CL upper limit on the branching fraction of 44 x 10~ 6 . 

We also recently measured the flavor-changing neutral-current decays D° — > e + e~ , D° — > 
and D° — > e ± {J. T that corresponds to an integrated luminosity of 468 fb~ l of data. To 
normalize the decays of D° — > £ + £~ , we use D° — > % + %~ control sample and applying a linear 
combination of Fisher discriminant [|]] of the following five variables: measured D° flight length, 
|cos0hei| angle between the momentum of the positively-charged D° daughter and the boost di- 
rection from the lab frame to the D° rest frame (all in the D° rest frame), the missing transverse 
momentum with respect to the beam axis, the ratio of the 2 nd and th Fox-Wolfram moments []|], 
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Dg->it + ^ + e" 



Figure 2: Invariant-mass distributions for Z) + — > 7t + £ + £^ (top) and Df — > 7C + £ + ^'' (bottom) candidates. 
The fit results shown in the solid line. 



and the D° momentum in the center of mass frame. To remove the continuum combinatoric back- 
ground we use the |co50hei| variable. Figure [7] shows distributions of |cr«0hei| before applying a 
minimum Fisher discriminant. 

The branching fraction of D —> £ + l~ is given by the following expressions: 



See ■ Nee 



where See is defined by 



and A^obs is denned by 



Iu/Su+Nbg- 



(3.1) 



(3.2) 



(3.3) 

7t + 7l~ candidates, respectively. The 



Nobs -- 

The Nee and are the number of D° -> and D° 
3$%% = (1.400±0.026) x 10~ 3 [10]. We use the likelihood ratio ordering principle of Feldman and 
Cousins [11] to determine 90% CL intervals. We find one event of D — > e + e~ with background 
of 1.0 ±0.5 events and two events of D° — > e ± /J. T with background of 1.4 ±0.3 events. These 
correspond to the 90% CL upper limits for the branching fractions < 1.7 x 10~ 7 for D° — > e + e~ 
and < 3.3 x 10~ 7 for D° — > e ± pL^. For the D° — > ;U + /i~ channel, we find eight events with expected 
background of 3.9 ±0.6. This corresponds to 90% CL upper limits on the branching fraction 
of [0.6,8.1] x 10~ 7 . Detail information on the likelihood selection, fitting procedure, systematic 



uncertainties, and fit results are available here [12 1. 



4. CONCLUSIONS 

We have searched for the decay modes Df\ — > % ± l T 6 

')+, D+ K+i+l^, and A+ ->■ 
pHfr £(!)+. No signals are observed and we report upper limits on 35 different branching ra- 
tios between 0.4 x 10~ 4 and 37 x 10~ 4 at 90% CL. This corresponds to limits on the branching 
fractions between 1 x 10~ 6 and 44 x 10~ 6 . 

We also have searched for the leptonic charm decays D° — > e + e~, D° — > pL + pL~ , and D° — > 
e ± fl T . We find no statistically significant excess over the expected background. These results 
supersede our previous results [ JT~3[ ] and are consistent with the results of the Belle experiment [|I4]]. 
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Figure 3: Invariant-mass distributions for D + — > K + £ + £^ (top) and Df — > K + £ + £^ (bottom) candidates. 
The fit results shown in the solid line. 



L D' >Ko'C j 


D + ^K"|i + |i + j 


D + -^K> + e + j 
^^^^ 


- D^K"e + e + 




Ds->K")i + e + 



1.9 2 1.8 1.9 2 1.8 1.9 2 

M(K eV) [GeV/c 2 ] M(K>V) [QeV/o z ] M(K"|i*s*) [GeV/c 2 ] 



Figure 4: Invariant-mass distributions for D + ^ K £ + £^ + (top) and D+ ->• (bottom) candidates. 

The fit results shown in the solid line. 
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Figure 5: Invariant-mass distributions for A+ — > p£ + &^ candidates. The fit results shown in the solid line. 

5. ACKNOWLEDGMENTS 

The author would like to thank the organizers of ICHEP 2012, the 36 f/i International Confer- 
ence on High Energy Physics, Melbourne, Australia. The supports from the BABAR Collaboration, 
the University of South Alabama, and the University of Mississippi are gratefully acknowledged. 

References 

[1] S. L. Glashow, J. Iliopoulos, and L. Maiani, Phys. Rev. D 2, 1285, (1970). 



5 



Charm Decays and Spectroscopy at BABAR 



Romulus Godang 



o 18 

? 1< 5 
S 14 

12 

r 10 

a 8 

in a 

1 \ 
o 



A c ->pe + e + 



A c ^pp. + p + 



A c ->pu. + e + 




2.25 2.3 2.35 

M(f5e + e + ) [GeV/c £ 



2.25 2.3 2.35 

M(pu + u+) [GeV/c^ 



2.25 2.3 2.35 

M(J5u+e + ) [GeV/c"; 



Figure 6: Invariant-mass distributions for — !> p£ + £^ + candidates. The fit results shown in the solid line. 
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Figure 7: Distributions of | cos^ei)! for the three signal channels: D° —> e + e~ (left), D° —> (center), 
and D° — > e ± jX zf (right). The bottom distributions show the signal Monte Carlo samples with an arbitrary 
normalization. 
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